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We propose a tomographic reconstruction scheme for spin states. The experimental setup, which 
is a modification of the Stern-Gerlach scheme, can be easily performed with currently available 
technology. The method is generalized to multi-particle states, analyzing the spin 1/2 case for 
indistinguishable particles. Some Monte Carlo numerical simulations are given to illustrate the 
technique. 
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I. INTRODUCTION 

The main idea of tomography is to reconstruct the 
density matrix g or, equivalently, the expectation value 
of any observable of the system from repeated measure- 
ments on an ensemble of identical states. In this pa- 
per a "spin tomography" for reconstructing spin states is 
developed in the framework of generalized tomography, 
starting from group theory iQ. There have been other 
proposals to infer the spin state . Our method is both 
easy to carry out experimentally and for the first time 
allows also a reconstruction of indistinguishable multi- 
particle spin 1/2 states, which is quite general since it 
encloses a vast class of experimentally accessible systems. 

The best known quantum tomographic procedure is 
optical homodyne tomography |^ , for the reconstruction 
of the density matrix gr of the radiation field from the ho- 
modyne probability p{x, 4>). It is based on the following 
formula ||] 



dx p{x, (j)) K{x — Xcf,) 



(1) 



where x^ is the quadrature operator and K{x) is an ap- 
propriate kernel function. We will not go into details 
about this formula, as we only want to stress the anal- 
ogy with the spin case. In fact, consider the spin density 
operator g, which is defined on a Hilbert space Tis of 
dimension 2 s -I- 1. We will prove the following formula 
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where the integral is performed over all directions of the 
versor n, p{n, m) is the probability of having outcome 
m measuring the self-adjoint operator s ■ n {s being the 
spin operator), and Ks{x) is a kernel function that will 
be defined later. 
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FIG. 1: Experimental apparatus for spin tomography. The 
Fizeau filter (F in the figure) selects particles with the same 
velocity from an incoming beam. These are then injected into 
a magnetic field Bi, which forms an angle ip with the y-axis 
and has intensity proportional to i9. For the tomographic re- 
construction the phases i? and ip need to be varied during the 
experiment. The remaining part of the apparatus is a conven- 
tional Stern-Gerlach scheme (in the figure we show the case 
of spin s = 1 as an example). A computer finally correlates 
the experimental results with the parameters i9 and p, in or- 
der to reconstruct the density matrix, according to formula 

(!)• 



It is possible to show that both Eq. (|^) and follow 
from a single operator identity, derived using group the- 
ory. In fact, define Tomographic Group G an unimodular 
group {i.e. left and right invariant measure are coinci- 
dent) which has a unitary irreducible square-integrable 
representation TZ{g) {g G G) on the Hilbert space H of 
the physical system. [A square-integrable representation 
is such that / dg \ {u\TZ{g)\v)\'^ < oo, where the integral 
is extended to all the elements g of the group G, dg is 
an invariant measure for G, and where the integral is not 
dependent on the choice of \v) € H (as will be shown 
later)] . The operator identity we derive (|l|) and (||) from, 
and which is valid for any tomographic group G, is the 
following 



dgTr[gTZ{g)] TZ^g) 



(3) 



valid for any trace class operator g acting on TL. This 
formula, derived from general considerations in ||l|, || is 
derived also in App. ^ using only group theory. 

The outHne of the paper follows. In Sect. ||, the to- 
mography procedure to reconstruct the spin state of sin- 
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gle particle systems is introduced and analyzed. The ex- 
perimental setup is described and some demonstrative 
numerical simulations of the procedure are studied. In 
Sect. [II the extension to the reconstruction of multi- 
particle spin states is studied. For distinguishable parti- 
cles, the general reconstruction procedure is given, while, 
for indistinguishable particles, the cases of two and three 
spin 1/2 particles is analyzed in detail. In Sect. IV the 
orders of magnitude of possible experimental setups are 
discussed showing the feasibility of the proposed method. 
In App. ^ the group tomography is derived by proving 
the tomographic reconstruction formula (^) in the frame- 
work of group theory. 



II. SINGLE PARTICLE SPIN TOMOGRAPHY 

Starting from the general operator identity (|^), we 
now specify the physical system as a single spin. In 
this case H — C'^''+^, s being the spin of the particle. 
For such a system, we can choose the group SU{2) of 
2x2 unitary matrices with unit determinant as tomo- 
graphic group G. In fact, SU{2) can be parametrized 
through the "rotation parameters" (n, ip) — where n — 
(cos sin I?, sin (/3 sin I?, cost?), £ [0,7r], £ [0, 27r], and 
ijj S [0, 27r] — and it induces a unitary irreducible rep- 
resentation on C^*+^. The operators constituting this 
representation are given by 



7^(r^, ^) 



(4) 



where s is the particle spin operator. Haar's invariant 
measure |^ for SU (2) is, with this parametrization and 
with the normalization needed for the invariant measure 
(see App. 



2s + 1 ip 
dg{n, ip) = ^ ^ sin^ — sini? dtp dd dip 



(5) 



As will be seen in the following, the choice of SU{2) as 
tomographic group G is not unique. It is easy to obtain 
the spin tomography (^ starting from Eq. (^, which 
now rewrites as 
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Evaluating the trace over the complete set of vectors 
\n,m) (which are the eigenstates of s • n, relative to the 
eigenvalue m), we find Eq. (0), by defining 
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and by noticing that (n, m|p|n, m) — p{n, rn). 

It should be pointed out that formula (m for optical 
homodyne tomography can be proven from Eq. (||) with 
a very similar argument. 



How do we use formula (g)? In order to measure the 
matrix elements gu — {ai\g\ai) for all i,l {{\ai)} being a 
basis for Tis), we only need to calculate (ai\Ks(m — s- 
n)\ai) and to measure p{n,m). 

The most convenient choice for the basis {|aj)} is the 
set {|r7i)} of eigenvectors of Sz {rn = — s, . . . , s). Thus, 
the calculation of the matrix elements of the kernel op- 
erator, by defining \i_m = m), yields 

{i\Ks{m-s-n)\l) (8) 
(2s -HI) 
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Observing that 
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with n±_ = (— sin c/j, cos (/?, 0), the evaluation of A/,m is 
given by 

Az^,„ = (;|e*''("'"'^*--™"'^*«V) = (^|e"''"^e-*''''''e*'^'*'-|TO) 
= e*'^("^'V(s + - m)\{s + iy.{s - ly. 



^ (-l)'^(cos |)2«+™-'-2^(_ sin D'-^+a-^ 
^—^ (s — Z — J^)!(s + TO — ^V-i'^ + ^ ^ to)!//! 



(10) 



where the sum is performed over the values of v for which 
the argument of the factorials is non-negative. In the last 
equality we used Wigner's formula [0. 



A. Experimental setup and state reconstruction 
procedure 

We now describe the method to measure the state g of 
an ensemble of non-charged particles, giving the details 
of the experimental apparatus, depicted in Fig. |l} The 
beam of particles impinges onto a Fizeau filter, which 
selects one velocity (in the x direction) for the particles. 
This is needed in order to ensure that each particle spends 
the same amount of time t in the subsequent region where 
a magnetic field Bi is present. The field Bi, which is 
parallel to the xy plane, is chosen so that Bi ~ Biti± ~ 
i3i(— sintp, costp, 0). In such way, its effect on the spin 
state g results in the unitary transformation g U , with 



U = exp [i7i3(sin (/3Sa; — cos (psy)t\ 



(U) 



Equation (11) follows from the Hamiltonian H = —jl- B, 
with fl = jhs {jl being the intrinsic magnetic moment of 
the particle, and 7 its giromagnetic factor). Successively, 
the particles cross a gradient of magnetic field B2, whose 
effect is to split the beam, giving a measure of Sz for 
the state Wg U, as in a, Stern-Gerlach experiment. In 
this way we obtain the probability {m\U^ g U\m), which 
is equal to p(n,m) by choosing Bi = — i3/(7i), and by 
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using Eq. (||). Therefore, by controlling the field Bi, we 
obtain p{n, m) for all n. In fact the direction of Bi selects 
Lp, while its intensity Bi selects "d. Now, in order to recon- 
struct the density matrix gi, only data analysis is needed, 
i.e. the insertion of the measured p{n,m) into Eq. (|^). 
One may object that an infinite number of measures are 
required. However, the calculation of the integral in (|^) 
with Monte Carlo techniques guarantees that the recon- 
structed matrix elements are affected by statistical errors 
only, which can be made arbitrarily small by increasing 
the number of measures. In practice, a rather small num- 
ber of data is required to obtain negligible errors, as we 
will show by numerically simulating the experiment. 
We first simulate the case of a coherent spin state [||, 




a e 



(12) 



FIG. 3: Diagonal elements of the matrix given in Fig. ^ The 
statistical error bars are obtained by dividing the measure- 
ments into 10 statistical blocks. The solid line indicates the 
theoretical value. 



where s± = Sj. ± isy. Notice the similarity with 
the customary optical coherent state, defined as \a) = 
goat-Q i|Q^^ where a is the annihilator operator for the 
optical mode and |0) is the vacuum state. In Figs. ||and^ 
we show the reconstructed density matrix g^oh = |a)ss (a| 
resulting from a Monte Carlo simulated experiment. 





FIG. 4: Density matrix for a thermal spin state gth- Here e = 
.75 and s = 2. A number of 60000 simulated measurements 
have been used in the reconstruction. 

B. Discrete spin tomography 



FIG. 2; Simulation of the reconstruction of the density matrix 
for a coherent spin state QcoH- The parameters for the state 
are a = 1 and s = 5. The simulation is performed using 3000 
spin measurements to generate the density matrix. 

As an additional example, in Figs. ^ and || we give the 
simulated reconstruction of a thermal spin state, which 
is the mixture defined by 



gth 



Tr[e- 



e e 



(13) 



The state gth describes a gas of non interacting spins in 
thermal equilibrium with a reservoir at a temperature T 
and in the presence of a magnetic field B^ parallel to the 
z-axis, i.e. e = —jhBz/{KBT), Kb being the Boltzmann 
constant. 



Up to now SU (2) has been used as tomographic group 
G for the reconstruction of the spin density matrix. This 
choice for G is not unique. For example, in the case of 
spin s = i, it is possible to use also the group defined as 
Q = {ia, —id, /, — /}, where d is the vector of Pauli ma- 
trices (T = {ax, (Jy,crz)- The following irreducible unitary 
representation on C ^ exists 



TZiiUa) — TZ{ — i(Ja) = CTa, a = X, fj, z 

n{i) = 7^(-/) = / . 



(14) 



Using this representation, from the tomographic recon- 
struction formula (BI) we obtain 



1/2 



(15) 



-1/2 a—x,y,z 



Notice that, by using Eq. (15) it is sufficient to measure 
the spin in only three directions. 
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FIG. 5: Main diagonal of the matrix given in Fig. The error 
bars, which in this case are practically negligible, are obtained 
by dividing the measurements into 10 statistical blocks. The 
solid line indicates the theoretical value. 



Analogously, for spin s = 1 it is possible to find a fi- 
nite group in alternative to SU{2). In fact, consider 
the 12 element tetrahedric group composed of the ±§7r 
rotations around the versors {ni — -^(1,1,1), n2 = 
^(l,-l,-~l),n3 = -i=(-l,l,-l),r?4 = 

of the TT rotations around {715 = (1,0,0), hq — 
(0, 1,0), fij = (0,0, 1)} and of the identity. It induces a 
unitary irreducible representation on the space C given 
by the 3x3 rotation matrices. Hence, Eq. (||) now be- 
comes 



7 

m— — 1 j — 1 



p{nj,m)ICj{m ~ s ■ fij) + -I , (16) 



with 



ICj{x) 



2 cos(27rx) 



J = 5,6,7 



(17) 



Notice that this procedure does not make use of a mini- 
mal set of measurements, since 14 experimental parame- 
ters must be determined in ([l6|), whereas there are only 
8 independent real parameters in the 3x3 density ma- 
trix. On the contrary, the case of spin s = ^ outlined 
previously does use the minimal set of measurements for 
such a system. In Fig. ^ a comparison between the two 
spin tomography procedures given by Eqs. (§) and (|l|) 
is shown through a Monte Carlo simulation. Notice that 
there is no significant difference in the results, showing 
that there is no substantial need for a procedure which 
involves a minimal set of measurements. 
For spins s > 1 an analogous procedure holds: one needs 
to find a finite group such that it induces an irreducible 
unitary representation onTi = C^^'^^. 



III. MANY PARTICLE SPIN TOMOGRAPHY 

The mathematical extension of the method to the case 
of a system composed of many spins is trivial, yet, it 
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FIG. 6: Monte Carlo comparison between continuous and dis- 
crete tomography. Continuous tomography uses SU{2) as 
tomographic group and is based on Eq. (p|), while discrete 
tomography uses 5(7(2) finite subgroups and is based on the 
reconstruction procedures given in Eq. ( p^ ) for s — ^ and 
Eq. (^|) for s = 1. Left: Convergence of the mean value of 
(sz) for a coherent a — 2 spin state for increasing number 
of experimental data (the theoretical value is given by the 
horizontal lines). The circles refer to continuous, the stars to 
discrete tomography. The upper graph is for spin s = |, the 
lower is for s = 1. Right: Plots of the statistical error bars of 
the graphs on the left vs experimental data. The error bars 
are obtained by dividing the experimental data into 20 sta- 
tistical blocks. Notice that the two tomographic procedures 
are essentially equivalent: they converge in the same way at 
the same result. 



predicts the necessity of performing measurements on 
single components and this may not always be possible 
when the system is composed of indistinguishable par- 
ticles. For this reason, we need to develop the theory 



A. Distinguishable spins. 



As tomographic group for a system of N spins we can 
simply use SU{2)^^ . Up to equivalences, its irreducible 
representations are given by the direct product of N op- 
erators (^ and the invariant measure is the product of 
N measures (H). As a consequence of the tomography re- 
construction formula (||) applied to S'C/(2)®^, we attain 
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readily the following generalization of Eq. 



(2.sfc + 1) '•2" 
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(18) 



where is the particle index. The trace term in ( [18|) gives 
rise to the probability p{ni,mi; • • • ; n^v, ttim) of obtain- 
ing TTik as result for the measurement of the fcth spin 
in the direction n^. This information is accessible only 
in the case of fully distinguishable spins. 

In Fig. ^ a simulated tomographic reconstruction of 
the value of (Sz) {Sz being the total spin component in 
the z-direction) is given for different multiparticle spin 
states. Notice how the number of the necessary experi- 
mental data increases exponentially with the number of 
spins, since the statistical error is exponential in the num- 
ber of particles. 
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FIG. 7: Left: Plot of (Sz) for different number of spins in a 
completely symmetrical state. A total of 10^ measurements 
for each mean value was performed in this simulation. Right: 
Semilog plot of the error bars vs the number of spins. Notice 
the exponential increase in the statistical errors. 



B. Indistinguishable spin 1/2 particles. 

Suppose we were given a system of N particles with 
the same spin. Such particles may be treated as identical 
by introducing a new dynamical variable, as in the case 
of the isospin. The spin density matrix (which is the 
partial trace over the orbital degrees of freedom of the 
global density matrix) is completely symmetrical, i.e. 



PgP- 



Q , 



(19) 



for any particle permutation P, because of the complete 
symmetry of the global density matrix. 

It is also possible to see that the spin density matrix is 
block diagonal in the representation of vectors of definite 
symmetry, with the subspace corresponding to each block 
spanned by vectors belonging to the same symmetry. In 
fact, given \4>) and 1-0) vectors of different symmetry type 



pL then {(plip) = 0. Hence, for any operator g, satisfying 
(|l9|), one has {(j)\g\ijj) = 0, as g\ip) belongs to the same 
symmetry type as {ip). 

Since the square of the total spin 5*^ and its z com- 
ponent Sz both commute with all permutation opera- 
tors P, the common eigenvectors of and Sz may be 
taken as a base for each of the diagonal blocks of the 
spin matrix. Let us now restrict our attention to s = 1/2 
spin particles. In this case, to each symmetry type there 
corresponds only one value of S, where S{S + 1) is the 
eigenvalue of S"^. In fact, given [A1A2] the partition of N 
which defines the class of permutations P that indicate 
a symmetry type, we find S — ^(Ai — A2) [ p^ . 

Let TLs^M be the space of vectors with assigned S and 
M [M being the eigenvalue oi Sz). The spin density ma- 
trix restricted to 'Hs,m, which is given by gs,M, is again 
completely symmetrical, hence [P, gs,M] — 0. Moreover, 
'Hs,M is associated with an irreducible representation of 
the permutations group |p^ . By using Schur's lemma, 
we can thus conclude that gs^M (x I, I being the identity 
in 'Hs,M- In T^s,m there may be vectors of different sym- 
metry type i, yet {i, S, M\gs^M\h S, M) does not depend 
on the index i, so that the probability for the measure- 
ment of S^ and Sz does not depend on the symmetry 
type. The same conclusion holds for the measurement of 
S^ and S -m for any versor rh. Hence, from the arbitrari- 
ness of TO, we conclude that blocks with the same S (and 
different symmetry type) are coincident. 

In conclusion, we have proved that in the {S*^} repre- 
sentation g is block diagonal, that each block corresponds 
to a value of S and that blocks with the same S are equal. 
Remarkably, applying Eq. (H) to each block, we can re- 
construct g measuring only the global quantities 5*^ and 
S ■ n. Some examples will clarify both the theory and the 
needed experimental setup. 

In the case of two spins 1/2, the spin density matrix 
will be of the form 



/ cm Cri2 CTis 

0'21 0-22 (723 

fSl <732 CTSS 

V a. 



>a , 



(20) 



where the a block corresponds to the subspace spanned 
by the eigenstates of S* = 1 (which are symmetrical with 
respect to particles permutations), while the a block to 
the subspace spanned by the only eigenvector of S' = 
(anti-symmetrical with respect to permutations). Ap- 
plying (0) to each block one finds 



dft 
in 



1,S -71 = M) K,=i{M - S -n) 



© piS^O). (21) 

According to (|2l|), in order to measure g, we only need 
the probability distributions p{S, 5 ■ n), corresponding to 
the operators S^ and S-n, for all n, which can be_suitably 
recovered using the apparatus depicted in Fig. 
will be analyzed later. 



which 
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Similarly, the spin density matrix of three spins 1 /2 is 
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(22) 



The ^ block corresponds to S* = 3/2, whereas the tt blocks 
both correspond to S = 1/2^ and are distinguished by 
their different symmetry properties. The argument pre- 



sented previously proves that 7^}^ - 
we can write = ^ © tt © tt, with tt 
applying (p|) to each block leads to 



vr?-, for all 



thus 



Again, 



e = 



Af=-f 

xK^^3{M - S -n) , 
xK^^i{M ~ S -n) , 



(23) 



(24) 



and the problem of determining g is again reconducted 
to the simultaneous measurement of S"^ and S ■ n. 

Both in the cases presented and in the general n spins 
case, the required experimental data are the distribu- 
tions p{S,S ■ ft). The apparatus to produce such data 
are basically equivalent in the two cases, as evident in 
Figs. ^ and ^, hence we shall limit the analysis to the 
two spins case. Here, the Fizeau filter and the magnetic 
field Bi = Bin± — _Bi(— sim^s, cos i^^, 0) have the same 
pvirpose as in single 
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FIG. 8: Experimental apparatus for the tomography of sys- 
tems composed of two spins i. 

Consider a beam of n non-interacting systems com- 
posed of two particles with spin 1/2. As the analysis can 
be immediately extended to a mixed case, for simplicity 
let us consider each system in the pure state 

= 7IO, 0) + 7^i|l, -1) + 7o |1, 0) + 7?|1, 1) , (25) 



with \a,b) standing for \S = a,M = h). The beam is 
split into three parts by the gradient B2 , and the systems 
arrive in detector B with a probability p{S = 1,M = 
1) = l7i P s-iid ill detector C with a probabihty p{S = 
1,M = —1) = |7"ip. The remaining particles reach 
position A with a probability 



and are left in the state 

= — a7lO,0)+7o"|l,0)). 

PA 



(26) 



(27) 



As the subsequent gradient is directed along the y axis, 
Eq. (|2^) is conveniently written using the eigenstates of 
Sy, i.e. \S,M)y: 

\^a) - — [7IO, 0), + 7o" -1), + 7o «o|l, 0), 

PA 



+7o"ai|l,l)y], 
where a. 



(28) 



y (l,i|l,0) {i — —1,0,1). Hence, the proba- 
tem to ar] 

[l7^P + KPl«on. (29) 



bility for a system to arrive at detector D is 
1 

PA 



PS 



By measuring and ps, the quantities and 
are obtained by inverting equations (|26| ) and (29). The 
coefficients |7o 17"^, |7i P and |7f:]^|^re the four prob- 
abilities p{S, M) we need for the reconstruction given by 
Eq. (|l|). 
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FIG. 9: Experimental apparatus for the tomographic recon- 
struction of the spin states of systems composed of three spin 



A similar argument shows that the equipment of Fig. 
^ supphes p{S, M), for all S, M, for a system constituted 
of three spins 1/2. 



IV. FEASIBILITY 

The orders of magnitude of the experimental param- 
eters are such that the experiment is feasible with cur- 
rently available technology. Only as an example, consider 
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the following cases of spin measurements of electrons or 
nucleons. For the magnet which is responsible for the 
field Bi with length of the order of 1 cm, we can mea- 
sure the state of a beam of electrons with speed ~ 10^ 
cm/sec, by using a magnetic field Bi = d/^t varying 
between and ^ 30 Gauss. On the other hand, in the 
nucleon case, choosing a speed of ~ 10^ cm/sec, we need 
Bi ranging between to 10^ 10'^ Gauss. Obviously, 
the parameters Bi and t can be adjusted over a wide 
range, according to the experimental situation. 



V. CONCLUSIONS 

We have presented a tomographic experimental pro- 
cedure for the measurement of the spin density matrix. 
The experimental scheme is a consequence of formula (|^) , 
which was proved using group theory. Through some 
Monte Carlo simulations, we have shown that the re- 
construction can be achieved with high precision using a 
limited number of measurements. The extension of the 
procedure to the reconstruction of states of multiple spin 
systems has been given, both for distinguishable spins 
and for indistinguishable spin 1/2 particles. Finally, we 
have shown that the orders of magnitude for the exper- 
imental setup are such that it can be implemented with 
currently available technology. 
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APPENDIX A: GROUP DERIVATION OF 
QUANTUM TOMOGRAPHY 

The proof of Eq. (j^) is obtained from the following 
lemma. 

Let A be an arbitrary trace-class operator on the 
Hilbert space TC of the system and TZ an irreducible uni- 
tary square integrable representation on TL of the tomo- 
graphic group G. Then 



= TZ{h) dg'n{g')\u){v\n\g'), 



TrA= / dgn{g)An^ig), 



(Al) 



where dg is an invariant measure for the group G, nor- 
malized as J dg \{u\TZ(g)\v)\'^ = 1 which is independent 
on the choice of the vectors \v) G Ti.. 
Proof: By using the unitarity of TZ and the properties of 
group representations, we can write for any h ^ G 



dgn{g)\u){v\TZ^{g)TZ{h) 
dg' Tl{hg')\u){v\n\g') 



(A2) 



which, through Schur's lemma, guarantees that 

dgn{g)\u){v\'R){g) = t„,„/„ , (A3) 



I-H being the identity in 71. Consider the quantity 

dg\{u\n{g)\v)\^ = T,,,, (A4) 



where \u) and \v) are arbitrary normalized vectors in 7i. 
From Eq. (A3) it is trivial to see that r^^y is independent 
on \u). One can check that it is also independent on \v) 
by noticing that, given an arbitrary vector \a) 

r,,, = I dgn{g)\v){v\TZHg)^ j dhn{h'^)\v){v\n\h-^) 

dhn{h)\a){a\n\h)^Ta,a , (A5) 

where the group unimodularity has been used in dg — dh 
with h = g^^ . Notice that the hypothesis of square- 
integrability of the representation guarantees the conver- 
gence of the integral in (A4). Thus, the natural choice 
for the normalization of the group's measure is to take 
Tu^t, = 1. The constant t^^v can be expressed in terms of 
■^11,1) by noticing that upon taking h = g^^ one has 



1=T, 

dh 



dg{b\n{g)\a){a\n\g)\b) = 
{a\'R{h)\u){v\n^{h)\a) v„ 



(A6) 



The lemma's thesis is now easily found by using the 
Schmidt decomposition oi A as A — ai\ui){vi\: 



dgn{g)AnHg) = Y.'^^ j dg n{g)\u,) {v,\n\g) = 

i 

Y,a^{v^\u^) ^TrA. (A7) 



Group Tomography Theorem. Let A be an arbitrary 
trace-class operator on the Hilbert space Ti of the system 
and TZ an irreducible unitary square integrable represen- 
tation on Ti of the tomographic group G. Then 



dgT,[ATZ{g)]TZ\g) 



(A8) 



Proof: Let O be an invertible trace-class operator, it fol- 
lows that TZ{g)0 is trace-class for any g E G. Hence it is 
possible to obtain, by applying (|A1|) twice 



dgTT[ATZig)]OTZHg)^ / dg' Tr[TZ\g')0]TZig')4A9) 
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Take a basis {|fc)} in 7i, one can obtain, using again 
lemma (Al), 



J dgTT[TZHg)0]{i\TZ{g)A\j) = 

f dg Y,{k\n\g)0\k){i\n{g)A\3) 
k 

Y,{k\OTr[\k){^\]A\J) = (^\OA\J) 



(AlO) 



From Eqs. ( [A9| ) and (AlO) it follows immediately that 
f do Tr[A7^(.g]]07^■^(g) = OA, which yields the thesis 
( [Aq ) by multiplying to the left both members by . 



It is trivial to extend theorem (A8) to the case of 
projective representations, i.e. group representations for 



which, given 51, 52, 33 G G such that gi ■ g2 = .93, one 
has 

7^(5l)7^(52) = e<(f-f^)7^(53) , (All) 

C G R being a phase factor depending on gi and (72- No- 
tice, moreover, that the theorem here presented is valid 
also for discrete and finite groups, with the sum on group 
elements replacing the integral. From result (A8), with 
an appropriate choice for the tomographic group and the 
irreducible representation, it is possible to prove the for- 
mula for spin tomography (H) -derived in the following 
section- and for optical homodyne tomography (|l|). No- 
tice that the unimodularity hypothesis given in the def- 
inition of tomographic group G can be relaxed without 
losing most of the results we give in this paper. 
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